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Abstract 

An effective theory approach is used to compute analyticahy the radiative cor- 
rections to the mass of the hght Higgs boson of the Minimal Supersymmetric 
Standard Model when there is a hierarchy in the masses of the stops (M^^ ^ 
Mj^ S> Mtop, with moderate stop mixing). The calculation includes up to two- 
loop leading and next-to-leading logarithmic corrections dependent on the QCD 
and top- Yukawa couplings, and is further completed by two- loop non-logarithmic 
corrections extracted from the effective potential. The results presented disagree 
already at two-loop-leading-log level with widely used findings of previous lit- 
erature. Our formulas can be used as the starting point for a full numerical 
resummation of logarithmic corrections to all loops, which would be mandatory 
if the hierarchy between the stop masses is large. 
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1 Introduction 



The Minimal Supersymmetric Standard Model (MSSM) predicts a light Higgs boson, 
with mass Mho of the order of the scale of electroweak symmetry breaking {Gp^ ~ f = 
246 GeV) times a small Higgs quartic-self-coupling. That Supersymmetry (SUSY) can 
naturally trigger this breaking, and stabilize the scale at which it takes place, is the 
most interesting part of the story (see for reviews and references). Here we take that 
for granted and our focus is on the perturbatively small coupling. Its smallness comes 
about because of two reasons: first, Supersymmetry dictates that the Higgs quartic self- 
couplings are given by gauge couplings (from D-terms) and by superpotential Yukawa 
couplings (from F-terms); second, the latter F-term contributions are absent in the 
MSSM since, in this model, quantum numbers prevent superpotential terms cubic in 
the Higgs fields. It is generic that quartic Higgs couphngs are directly related to 
the Higgs mass after electroweak symmetry breaking (the Standard Model is the best 
known example). All this results in the well known tree- level upper bound M^o < 
M| cos^ 2/9 [the SUSY parameter tan/5 is the ratio V2/V1 of the two Higgs vacuum 
expectation values (vevs). We follow the usual convention: V2 generates the mass of 
the top quark and Vi that of the bottom quark]. 

Radiative corrections to M^o can be quite important because of top-stop loops that 
introduce a dependence on the top Yukawa coupling, ht, which is sizeable, while this 
coupling does not enter in the tree-level Higgs mass. This can lead to cases in which 
one-loop radiative corrections to Mf^o are comparable to, or even larger than the tree- 
level part of it (without this being an indication of the failure of the perturbative 
expansion) . 

In addition, the one-loop corrections to M^o are logarithmically sensitive to the mass 
ratio, mi/mt, of the average stop mass over the top mass, which could be large if there 
is a hierarchy, Msusy / ^ew ^ 1, of the SUSY mass scale over the electroweak scale. 
As a consequence, radiative corrections to Myfi beyond one-loop can be important if 
hi{mi/mt) is large. In that event, standard renormalization group (RG) techniques can 
be used with advantage to resum these logarithmic corrections to all loops. 

During the last decade, the precise determination of M^o as a function of the super- 
symmetric parameters has received continued attention The development of 
increasingly refined calculations of Myfi is the story of a stepwise climbing of this ladder 
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of loop corrections and has been told elsewhere (see e.g. for a brief account) so it 
will not be repeated here. The current status of what has been achieved, by the com- 
bined use of direct diagrammatic calculations, effective potential methods and RG tech- 
niques, could be summarized in this way: all one-loop corrections are known and 
the dominant two-loop corrections of order (9(asQ;tm^) and 0{a'^m^) are also known, 
including finite (non-logarithmic) contributions [here as = (7^/(47r) and at = /i^/(47r), 
with Qs the QCD gauge coupling]. Higher order corrections at leading-log and next- 
to-leading-log order [O^atmla^J^n^mi/rrit)]^) and 0(atmja"^^[ln(m(-/mi)]")] can be 
resummed using one-loop and two-loop RG /^-functions, respectively. 

Beyond tree level, Myfi is sensitive to many SUSY parameters, but the most impor- 
tant are those of the stop sector (and of the sbottom sector also for large tan/5). They 
are given by the stop mass matrix: 



M 



Ml + ml mtXt 
nitXr Ml + 



(1) 



where we have neglected D-terms, M| (M|,) is the soft-mass for ii (iji), and 

= A + /i7tan/3 , (2) 

with At the soft trilinear coupling associated to the top Yukawa coupling and fj, the 
supersymmetric Higgs mass in the superpotential. 

The dependence of the radiative corrections to M^o on these parameters has been 
studied before in different specific regimes. In this paper, we focus on the case in which 
there is a double hierarchy. Ml ^ M^j ^ rrit (the case Ml ^ Mr ~ rrit can be worked 
out along similar lines). In this situation one should care, not only about potentially 
large logarithms like the usual ln.{ML/mt) and ln(Mij/mt), but also about ln(M2,/Mij). 
Radiative corrections to M^o for this type of stop spectrum have been considered in the 



past [|ll| , pr2| but there is room for improvement, as we will show. First, if the hierarchy 
between the stop masses is large, a numerical resummation of logarithmic corrections to 
all loops is necessary to get an accurate determination of the Higgs mass and, in order 
to do this, one has to identify first the relevant RG functions and threshold corrections. 
So far, this has not been done. Second, although previous analyses represent important 
steps ahead, they are not complete in one sense or another: either they do not include 
all potentially relevant corrections or, if they do, the corrections are not cast in a form 
suitable for RG resummation. 
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The plan of the paper is the following. In section 2 we present the main calculation. 
We use an effective theory method to extract and classify all two-loop dominant {i.e. 
ht and (^s-dependent) radiative corrections to M^o. The result of this calculation can 
be used as the starting point for a full numerical evaluation of the Higgs mass in the 
case of hierarchical stop spectra, although we do not undertake that task in this paper. 
In section 3, we discuss the possibility of finding a one-loop 'improved' approximation 
to Mho that, playing with a judicious choice of the scales at which parameters are 
evaluated, tries to absorb higher order corrections. This exercise is a good point at 
which to compare our main result, presented in section 2, to previous analyses existing 
in the literature, with some of which we disagree already at the level of two-loop 
leading-log corrections. We dedicate section 4 to such comparisons. Section 5 presents 
our conclusions and outlook for future work. For reference. Appendix A presents 
an explicit formula for M^o which includes up to two-loop-next-to-leading logarithmic 
corrections. Appendix B is devoted to the calculation of two-loop threshold corrections 
for the Higgs quartic self-coupling, of direct interest for the completeness of the two-loop 
calculation of M/jO. Finally, Appendix C gives the relationships between MS running 
parameters (in which our results are expressed) and on-shell (OS) quantities. 

2 Effective theory calculation 

We consider the MSSM with a particle spectrum in which all super symmetric particles 
have a common mass, Msusy, much larger than the electroweak scale (say a few TeV) 
except for the lightest stop, which is much lighter although still heavier than the top 
quark. In particular, we remark that the mass of the pseudoscalar Higss, M^o, is also 
taken to be Msusy, and therefore, the model contains just one light Higgs doublet. To 
be precise, and referring to the stop mass matrix written in eq. (|T]), we consider 

mt < M« < Mi = Msusy ■ (3) 

Concerning stop mixing, we also assume that it is not too large, so that it is a good 

approximation to say that the lightest stop is mostlyQ t^, while the heavier one is 

mainly tu- In other words, we are in a situation in which the stop mixing angle is small. 

Nevertheless we do keep the dependence with the stop mixing parameter Xf and we 

^ This avoids problems with a large contribution to Ap, which would be present in the opposite 
limit in which the light stop is mostly t]^. 
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will derive our results as a series in powers of mtXt/Ml (note that our approximation 
is nitXt/Ml <^ 1, not Xi/Ml ^1). The case of a hierarchy in stop masses due to very 
large Xt (rather than to different diagonal soft masses) is worth separate study but it 
is more complicated and we do not consider it here. 

To compute the radiatively corrected Higgs mass in the hierarchical case (|^), we 
make use of an effective lagrangian approach, descending in energy from Msusy down 
to the electroweak scale mj. In doing so we encounter different effective theories at 
different energy scales. Above Ml = Msusy the relevant theory is the full MSSM. 
Between Ml = Msusy and Mr the effective theory contains only the Standard Model 
particles with a single Higgs doublet (that particular rotation of the two Higgs doublets 
of the MSSM which is responsible for electroweak symmetry breaking and has SM 
properties) and, in addition, the light stop. Below the mass scale Mr of that light 
stop the effective theory is simply the pure Standard Model (with calculable non- 
renormalizable operators, remnant of the decoupling of heavy SUSY particles). 

To compute M^o we start at Msusy with the known value of the quartic Higgs 
coupling. A//, as a boundary condition fixed by Supersymmetry. We run this coupling 
down to rrit in the different effective theories just mentioned, taking care of threshold 
corrections whenever some energy threshold is crossed. The procedure is standard and 
follows the general prescriptions for effective theory calculations. For general reviews of 



this subject we refer to and references therein. We also found useful some general 



discussions in ref. ||T9[, a more specialized paper which studies the effective theory of 
a linear 0{N) sigma model. Similar effective theory techniques have been applied to 



study the decoupling limit of the MSSM with heavy superpartners [20 



We work in an approximation that neglects in radiative corrections all couplings 
except Qs and ht [our results could be extended easily to include also hf, (bottom- 
Yukawa) corrections, which can be significant for large values of tan/3]|. We keep 
electroweak gauge couplings only in the tree- level contribution^ to M^o. In this con- 
nection, the quartic Higgs coupling, X^, is considered to be itself of one- loop order 
[Xh ~ ^t/(167r^)] when it appears in radiative corrections. Within this approximation, 

•^Following this approximation, we neglect the effects of gauge couphngs in the masses of SUSY 
particles (which only affect A/^o through radiative corrections). In particular, Higgsinos simply have 
mass \n\ = Msusy- 

•^For numerical applications the full dependence on gauge couplings at one-loop is known and can 
be included. 
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we plan to extract analytically the radiative corrections to M/jO up to two-loop order, 
that is, we compute one-loop leading-log and finite terms plus two-loop leading-log, 
next-to-leading-log and finite corrections to M^o . We also use, whenever necessary, ex- 
pansions in powers of the mass ratios nit/ Mr, nit/ Ml and Mr/ Ml, as is common use in 
effective theory calculations. If these ratios are not small there is no necessity of using 
effective theory methods: the corresponding corrections to M^o are not logarithmically 
enhanced and other existing calculations should be valid. 

The analytical result for M^o that we obtain in this way is interesting for two 
reasons: first, it identifies the ingredients (threshold corrections and renormalization 
group functions) necessary for a full numerical computation of the running of Xh'-, 
second, it is interesting in order to develop simple and compact approximations to 
the full numerical results. Our goal is then to find a two-loop formula for the Higgs 
mass in which the corrections are classified in such a way as to permit a numerical 
resummation of leading and next-to-leading logarithmic corrections to all loops. As 
explained, this resummation is mandatory if the hierarchy (^) is sizeable, when simple 
analytical approximations start to fail. We defer that numerical evaluation of M^o to 
a future publication and concentrate here upon the analytical study. 

2.1 Plan 

The method we follow is very similar to that used by Haber and Hempfiing in ref. 
to compute radiative corrections to M^o for low values of the pseudoscalar mass, M^o, 
case in which the theory below Msusy is a two-Higgs-doublet model. The plan of our 
calculation is to integrate the equation dXu / dhiQ"^ = [3\^ from Msusy = Ml, where 
Xh is related to SUSY parameters, to Mew = ^t, where determines the Higgs mass. 
Taking into account the different running in the two effective theories, above and below 
the intermediate threshold at Mr, and writing explicitly the threshold corrections to 
Xh [with a superindex, as in X%{Ml), we always indicate a tree- level value], we find 

X{mt) = XUMl) + SXh{Ml) - /3x{Q)dlnQ^ 

jQ=mt 

+5XH{Mn) - /3,^(Q)rflnQ2 . (4) 

JQ=M+ 

The quantities 5Xh are the threshold corrections for Xh at the indicated scales. We 
call A the Higgs quartic coupling below M/j to distinguish it from Xh above Mr. 
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If we next expand the /3-functions around a particular value of the scale, and make 
a loop expansion up to two-loops [(3 = P^^^ + P^"^^ + ...] we get 



Kmt) = XHiML) + 5MML) + S\HiMn) 



2(ilnQ2 ml 2rflnQ2 Ml 



rrit 



In 



Ml 



It is important to make explicit the scale at which one-loop /3-functions are evaluated, 
because different scale choices amount to a two-loop difference [the scale choice in two- 
loop terms like /J^^) or djS'^^^ / d\nQ'^ has only effects starting at three loops]. The same 
comments apply to the choice of the scale at which to evaluate the masses inside one- 
loop logarithms. In (H), the RG procedure dictates that they are evaluated at a scale 
equal to the mass itself, that is, M/j = Mji{Mji), = Ml{Ml) and rrit = mt{mt). In 
a similar way, it is important that the couplings which appear in /?*^^)'s are evaluated 
taking into account the corresponding one-loop threshold corrections. All this will be 
shown more explicitly in the following subsections. 

Eq. (^) already illustrates some properties of radiative corrections which are generic: 
i) Leading- log contributions at any order depend only on one- loop RG-functions and 
are, therefore, insensitive to threshold corrections and two-loop or higher RG-functions. 
The reason is simple: by definition, in leading-log corrections each power of the loop 
expansion parameter a [at or in our case) is accompanied by a logarithm (that 
arises from RG running between two mass scales). However, threshold corrections in- 
troduce powers of a's without such logarithms, while n*^-order RG-functions introduce 
a factor a""^ for each a log. ii) Next-to-leading-log terms are instead sensitive first to 
two-loop RG-functions and second, to one-loop RG-functions times one-loop threshold 
corrections. In turn, they are not sensitive to three-loop (or higher) RG-functions or 
two-loop (or higher) threshold corrections. 

With this hierarchical classification of radiative corrections in mind, our calculation 
aims at finding the relevant one and two-loop RG-functions plus one-loop threshold cor- 
rections. This would allow the resummation of leading and next-to leading logarithmic 
contributions to M/jO to all loops. Nevertheless, in our analytical formulas we stop at 
two-loops, including also two-loop non-logarithmic terms. 
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2.2 SUSY threshold: matching MSSM with SM + 

The effective theory below the supersymmetric threshold at Msusy is described by the 
most general Lagrangian built with SM particles plus tn, non-renormalizable in general 
but invariant under the SU{3)c x SU{2)l x f/(l)y gauge symmetry: 

^SM-,t, = {V,inrHvnnr + iV,h)\V'^h)-Ml\in^^ 



_ 2 2 

->^'Huii%*h) . d^Ml) - K'W\h\'' - x\h\'\iR\' + ... , (6) 



where the ellipsis stands for terms of higher order both in fields (bosonic or fermionic) 
and derivatives. The Higgs doublet field is represented by /i, is the top-bottom 
quark doublet and the right-handed top quark field (a is a colour index). The dot 
(.) stands for the SU{2) invariant product and c = — zo"^. We have written explictly 
only the third generation Yukawa coupling, which in this intermediate-energy theory 
we call Qt- We keep only terms directly related to our calculation and do not write, for 
example, fermion kinetic terms. 

The parameters in the Lagrangian (^) are determined by matching at the scale 
Q = Msusy with the full MSSM theory, i.e. by requiring that the effective theory 



and the full MSSM give the same physics at low momentum ||T8[. To do this matching 
at tree level, we first obtain the equations of motion of the heavy MSSM fields and 
substitute them (making a low momentum expansion) in the MSSM Lagrangian. What 
one obtains is (we use a prime to distinguish this Lagrangian from that of the MSSM 
with heavy fields not replaced by their equations of motion) 



2 



^'mssm = (2^;.tR)"^(2^^tR)" + (2^^/^)^(2^'^/^)-M^|t^|2-rn2|/i| 



-T2' 



_ui„4:\Y |2 



iRh-PL{d')h h.P,{d')hil 



htclslX,Y:{il*h) . PL{d'rnPA{d')h\in\' + h.c. 

+ htclsl\X,Y,nr^*h) . PL{drtRPA{drt'RPL{d')'hfR 
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-hlclPA{d'){Q"^,ctr){QlJc tf) + ... (7) 

where g and g' are the SU{2)l and U{1)y couphng constants respectively; h is the 
SU{2) conjugate of the hght Higgs field h; cp = cos/5, S/3 = sin/3; = (1, —a); T^^ 
are the SU{3)c generators in the fundamental representation; and 

Yt = At- ^* tan /3 . (8) 

We have written and in Csm+Ir distinguish them from M|, and rn? in 
^'mssm- W^e have also introduced the operators 
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Ml + 92 



(9) 



for X = L,A,H,G (corresponding to the soft mass of t^, the pseudoscalar mass 
Ma, the higgsino mass = and the gluino mass Mq, respectively). In eq. 
the subindex E indicates a low-momentum expansion in powers of d'^/M^. In (|^), these 
operators act only inside the square brackets they are in. 

The origin of each non-renormalizable term in C^jg^j^i [eq. (|^)] is easy to interpret 
as coming from the tree-level exchange of one or more heavy particles [identified by 
the propagator operators Pxid"^)]- This is shown diagrammatically in figure |I], which 
shows the tree-level diagrams that give rise to the different terms of (|^, upon col- 
lapsing heavy particle lines to a point. The line code we use is the following: a thin 
dashed line with a small arrow [which indicates the flow of SU(2) quantum numbers] 
represents the light Higgs doublet; the same type of line but with double dash cor- 
responds to the heavy Higgs doublet; a double continuous line with an SU{2) arrow 
represents a Higgsino; a dashed bold line with a large arrow [indicating SU (3) colour 
flow] represents a light stop; the same type of line but thicker and with a larger arrow 
[which indicates the flow of SU (2) and SU (3) quantum numbers] is used for the heavy 
stop; gluinos are represented by a continuous line with a wiggle; a top-bottom quark 
doublet is represented by a solid line with a large arrow, while the same type of line 
with a smaller arrow corresponds to a right-handed top quark. For our calculations 
and diagrams we work in the unbroken-symmetry phase, with the full SU{2) doublet 
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x^ X, ^ 




Figure 1: Tree-level diagrams, with heavy supersymmetric particles interchanged, that give rise to 
different terms of the Lagrangian (|^) . 

structure unresolved. This simplifies our task: when dealing with separate diagrams for 
the different components of a doublet there are cancellations, due to SU{2) symmetry, 
which are immediately obvious in our approach. One such example is the diagram of 
figure 1^. There is a cancellation between the contributions of different SU (2) doublet 
components running in the loop. If one works with complete SU{2) multiplets, this 
cancellation is reflected in the impossibility of drawing such a diagram properly: there 
is no choice for the arrow of the heavy Higgs field in the loop that is consistent with 
the flow of SU (2) quantum numbers through other lines of the diagram. 

Comparing C\,[ggj^,^ [eq. (|^)] to /^sA/+tj{ [^Q- ®\i we get the tree-level matching 
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Figure 2: One-loop diagram that cancels when all the components of SU{2) doublets run in the 
loop. 



conditions: 







Xhu (Ml) 


P^t P 1 2 

= 9t 9t cos2(3 + 6Xhu 


MMl) 


= lg! + ^g" + 6Xu{ML) , 




= htsm(3 + 6gt{ML) , 


mI{Ml) 





(10) 



where we have just indicated the presence of loop corrections and, with an abuse of 
notation, used gt = ht sin f3 everywhere. In addition, we hst the following tree-level 
threshold values for some non-renormalizable couplings in Cgj^f^i^^ which will play a 
role in higher loop calculations: 



9t^ + SX'Hu{ML) 



9t- 



Mi 
Mi 



Ml 



+ 



Ml 



1 - 



XtY: 



Ml 



'IV. 



This tree-level matching would be enough if we were only after leading-log correc- 
tions to M/jO (which are not sensitive to threshold corrections). As we have discussed 
already, if we want to correctly obtain next-to-leading log contributions, this matching 
must be done at one-loop level. That is, we need the one-loop threshold corrections 
in the matching conditions (PH]). To compute them we match the one-loop effective 
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Figure 3: One-loop diagrams for the threshold correction to Xh in the full MSSM. 



actions in both theories, MSSM and SM + t/j. Again, to do this, we first substitute in 
the MSSM effective action the equations of motion of heavy fields in a local momen- 
tum expansion. In other words, we match ILPI (one-light-particle-irreducible) graphs 
with light-particle external legs in both theories. This procedure leads to the one loop 
threshold corrections: 



SoXh{Ml) 
SoXhu{Ml) 



SogtiML) 

Here, N, 



87r2 
87r2 



c 4 

9t- 



Mt 



1 + + O 
Ml 



'M£ 



1 + 2Mg 



Xt + h.c. 



Ml 



4 + (1 + 2si + N, 



/3 



Ml 



+ s 



Y,x: + h.c. 

Ml + 



AMI 



3c 



^ Mt 

\XtYt\'' 
/3 4M| 



O 



'M^ 
Ml 



''''-C,{Ny-^X, + 



Mr 



9. 



87r2 

2 



Ml 



'Ml 
Ml, 



9t 



(12) 



(13) 



(14) 
(15) 



47r2 * 87r2 

3 is the number of colours and C2{Nc) — {N^ — 1)/{2N^ is the quadratic 
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Figure 4: One-loop diagrams for the threshold correction to the quartic Higgs coupling Xh in the 
effective theory SM + tj^. 

Casimir of the fundamental representation of SU{Nc). For simplicity in our notation, 
we remove the bar from in these equations and from now on, in the understanding 
that Mji is, in what follows, the parameter of the intermediate-energy theory. 

These results are approximations in which we have neglected all couplings other 
than the top Yukawa coupling and the strong gauge coupling. They are expansions 
in powers of Mr/ Ml to the indicated order of approximation. We have truncated the 
expansions in such a way as to reproduce the correct one-loop corrections to M^o up 
to and including 0{M'j^/ Mf^) terms while we neglect this type of corrections in the 
two-loop contributions [for that reason we need to keep (9(M^/M|) corrections only 
to \hi which enters M^o already at tree level]. In addition, we have used the assumed 
identity of different heavy masses {Ma = Mfj = Mq = Ml) to simplify the expressions 
(although we leave explicit linear terms in Mq to allow for possible sign effects related 
to that mass). 

In figure ^ we give the ILPI diagrams that contribute to at one loop in the 
full MSSM, while figure § shows the corresponding diagrams in the SM + theory. 
Couplings in this last figure are distinguished from those in figure ^ by a black square 
to represent that they already include tree-level matching corrections. The one-loop 
threshold correcction 5\h is given by the contribution of the diagrams of figure ^ minus 
the contribution of the diagrams of figure ^. We therefore omit from these figures those 
diagrams that would be exactly equal in both theories (such diagrams do not contribute 
to the threshold correcction 5\h) or diagrams that are simply zero. 
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Figure 5: One-loop diagrams for the threshold correction to the Higgs-stop coupling Xhu in the full 
MSSM. 
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Figure 6: One-loop diagrams for the threshold correction to the Higgs-stop coupling Xhu in the full 
MSSM (continued). 
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Figure 7: One-loop diagrams for the threshold correction to the Higgs-stop coupling Xhu in the 
effective theory SM + in. 

In a similar way, figures § and ^ give the ILPI diagrams that contribute to Xhu 
at one loop in the full MSSM, while figure |^ shows the corresponding ILPI diagrams 
in the SM + t/j theory. At the order we work, the only diagrams that contribute to 
the threshold corrections of the top Yukawa coupling, gt, and of the mass Mf> of the 
light stop tn, are diagrams in the full MSSM. They are shown in figure |]for gt, and in 
figure ^ for Mr. 

The subindex in eq. (|T5D is meant to indicate that these are not the final one- 
loop threshold corrections: there are also threshold corrections to kinetic terms and, 
after redefining the fields to get canonical kinetic terms, the results in eq. ( [T^ ) are also 
affected and one finally gets 



Si\h{Ml) 
6iXhu{Ml) 



6oXhu{Ml)-Xhu{SZu + SZh) , 
6ogt{ML) - ^gt {SZ^^ + SZ^^ + 6Zh) + 
5oMl{Ml) - Ml5Zu , 



(16) 



with wave-function threshold corrections encoded in 



5Zi 



647r2 



(3 + 4) , 



15 




Figure 8: One- loop diagram for the threshold correction to the top Yukawa coupling in the full 
MSSM. 




Figure 9: One-loop diagrams for the threshold correction to the mass and wave- function of Ir in 
the full MSSM. 




Figure 10: One- loop diagrams for the threshold correction to the wave-function of the light Higgs 
h in the full MSSM and the effective theory SM + tR. 
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Figure 11: One- loop diagrams for the threshold correction to the wave- function of the left-handed 
top quark in the full MSSM and the effective theory SM + tR. 





Figure 12: One-loop diagrams for the threshold correction to the wave- function of the right-handed 
top quark in the full MSSM and the effective theory SM + tR. 



^9t I V 12 



,2 



1 



327r2 

SZ, . _|i,C,(iV.) + ^ (l + + . (17) 



for Q1, t% h and fields, respectively. We have also added a term (see pll ' 



9 9^ 

^DR^MS = 16^'-^2(A^c 



to correct from the change of scheme, from DR ' (the modified DR scheme of [|^) in 
the MSSM to MS , which is the scheme we use below the supersymmetric threshold. 
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Figure 13: One-loop diagram for the threshold correction to the top Yukawa coupling in the full 
MSSM with light-heavy Higgs doublet mixing. 

The diagrams that contribute to the threshold corrections to kinetic terms, given in 
(p!7|), are shown in figure ^, for t^; in figure IC for h; in figure ^ for ti and in figure 12 



for tji. In the last three figures we give together the diagrams in the full theory and 
(with a minus sign in front) those in the effective theory (SM -|- t/j). 

Several comments on the one-loop threshold corrections we have presented are in 
order. To get the correct matching conditions it is important that both heavy and 
light particles propagate in loops when computing the full MSSM effective action (see 
e.g. the discussion in |T^). We have computed the above threshold corrections evalu- 
ating the functional determinant expression for the effective action and also by direct 
diagrammatic calculation of the matched graphs (for this task, the general reference 
23| was helpful, as usual). The former method has the advantage of being systematic 



and of simplifying the determination of symmetry factors, the second illuminates the 
physical origin of what is being computed. We find agreement between the results from 
both methods. 



As expected on general grounds |T8[, the threshold corrections are infrared finite, 
i.e. all sigularities in the limit Mji cancel in the matching. This happens for 
the threshold corrections and for all their derivatives with respect to the light mass 
Mr. In particular, no dependence on InM/j is left in threshold corrections. As is well 
known, for this cancellation of infrared divergences to occur, it is crucial to keep enough 
derivatives in the low-energy effective couplings of CgMj^f^ [eq. (|)]. This successful 
cancellation provides a partial check of our results. 

The treatment of the matching in the Higgs sector requires some explanation. For 
large pseudoscalar mass, M^o, we can rotate the original two Higgs doublets of the 
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MSSM, and H^, into light and heavy doublets {h and H, respectively) which are 



given by 

h \ ( cos [3 —sin 13 




(19) 



H j \ sin /5 cos B 

At one- loop order, however, a mixing between h and H is induced {e.g. by stop loops). 
As we would like our field h to be the true light Higgs doublet at one-loop, we treat the 
Higgs sector at one-loop in the matching of effective theories. More precisely, we use 
one-loop equations of motion for the heavy Higgs field. If this is done, all ILPI diagrams 
in the full theory that contain a external leg in which H switches to h through a loop 
correction, cancel exactly with the one-loop contribution from the equation of motion 
for H. This is what one expects of Higgs fields properly diagonalized at one-loop level. 
As a result, such diagrams do not contribute to one- loop threshold corrections and, 
for that reason, we have not included them in previous figures. An example of such 



diagrams is given in figure [T^ for the top Yukawa coupling. 



2.3 Running down to Mji 

Once we have fixed the couplings of the intermediate-energy theory at the scale Q = 
MsusY in terms of the parameters of the full theory, we run them down in energy until 
we reach the next threshold at M^. The renormalization group equations (RGEs) in 
the intermediate effective theory described by Cgj^^^i^ [eq. (|)] can be easily computed 
at one-loop {e.g. through the effective action) including where necessary the effect of 
non-renormalizable operators. For two-loop results we particularize to this theory the 



general formulas presented in [24 



Following the notation introduced in eq. (^ for the couplings of the intermediate- 
energy theory, we have, for the Higgs quartic coupling: 



167r2 



fl67r2)2 



\gl^1glC,{N:){\\^-g\)-\\^ 



(20) 



with 



■IS = -«^iL . (21) 



IGtt 

describing the wave-function renormalization of the Higgs field 
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In view of eq. we need to compute d/^l^/rf In and, therefore, we also need 
the RGEs for Xhu and gt (albeit only at one-loop). Those for and X'^jj will not be 
needed because of the additional factor M|, which we neglect in two-loop corrections. 
For the Higgs-stop coupling we find 

/^S, ^ T^^^^ [(^'^ + 1)^^ + 2^^^ + 3^^] - ^HU^^H + lu^) , (22) 

with 

7j;'^ = (3-a)-|i^C2(iV,), (23) 

describing the wave-function renormalization of the stop field (in Landau gauge, in 
which we work, a = 0). For the top- Yukawa coupling we get 



0(1^ ^("9 2_8„2, .24) 



Finally, we also give 



^^i"^ilnQ2- evr^^^' ^^^^ 



which will be useful later on to make contact with effective potential results. We remind 
the reader that Mii{Q) is the running mass of Ir in the intermediate effective theory, 
not in the full MSSM. 

These RGEs are approximations in which the SU{2)l and U{1)y gauge couplings 
and all Yukawa couplings other than gt are neglected. We keep Xh in (^) because it is 
formally of order gf/{lQ7i'^) and, therefore, the XhIh^ term in (3^^^ is of the same order 
of f3^^ and should be keptQ. 

In we also keep the contributions from non-renormalizable operators (k^ and 
^'hu) "which are suppressed by a M'^/Ml factor. This is necessary if we want that 
precision in one loop corrections to M^o. At two- loops we do not keep such subdominant 
terms and, consequently, such contributions have not been kept in or Px^^- 

2.4 Intermediate threshold: matching SM + Ir with SM 

With the RGEs presented in the previous section we can then run the parameters down 
to M/j. At Mji we need to match the low-energy effective theory, which is the SM (with 

^In fact, XhJh"^ is even dominant with respect to /3^^ because of the logarithmic enhancement 
Xh ^ gf /{IGtt'^) ln{ All/ mf). This term must be included even if we were only interested in leading- 
log corrections . 
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non-renormalizable operators formed out of SM fields): 

O 1 , lO 1 



A|/Lr + M2ct5?-A + h.c.] + ... 



(26) 



to the intermediate effective theory, SM + t/j, given by the Lagrangian (|]) but with the 
heavy field tji removed. Due to the fact that no low-energy SM fields (or combinations 
of them) have the quantum numbers of t^, there are no threshold corrections at tree 
level to the couplings in (pBI) . At one-loop we find no threshold correction for the top 
Yukawa coupling, 5ft, and a non-zero correction to the Higgs quartic coupling: 



X{Mr)=Xh{Mr) + 5X{Mr) 



(27) 



with 



SX{Mr) 



87r2 



Nr 



87r2 



9t 



k\mI) - Xhu{MI)X'h^{MI) 



Ml + O 



Mi Ml 



+ 



(28) 



The diagrams that contribute to this correction from the intermediate theory are those 
already depicted in figure |. In (0), the couplings that appear in the first expression 
should in principle be evaluated at the scale Mr. At the order we work, however, the 
choice of that scale is not important and we can simply use their tree-level values as 
computed at Ms and given in eqs. (|10|) and ( pJ]) . In this way we get the final expression 
in terms of SUSY parameters. The difference between this expression and the proper 
one is a two-loop next-to- leading- log term of order M|./M| and we neglect such terms. 



2.5 Running down to nit 

Next we use the SM RGEs to run A and gt from Mr down to the electroweak scale, 
say to Q = rrit. These RGEs are the following. For the Higgs quartic coupling, A 
[normalized to M^o = Xv"^, with v = 246 GeV] we have 



167r2 
(167r2)2 



(29) 



with 'j^j^^ as given in eq. (PTf). We again keep the term X'y^jp in order to get the 
correct two-loop result for M^q, but neglect all couplings other than gt or gg. For 
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the top Yukawa coupling we get the same one-loop running as in the intermediate 
theory, eq. (^^. Notice that the RGEs for the intermediate theory, eqs. (PO]), reduce 
to eqs. in the limit k^, Xhu ^ 0, as they should. 



2.6 Mass formula 



Once we have \{mt) we simply use the SM relation 



1 + 



9t 



to obtain the Higgs mass. Here v is the Higgs vev, running with RGE 



(30) 



(31) 



rflnQ2 

with as given, at one loop, in (pTj). The one-loop correction factor in (pO]) takes 
care of Higgs wave-function renormalization effects I^Sj. The fact that A itself is to 
be considered of one-loop order makes it unnecessary to refine eq. ( pO]) with two-loop 
effects. 

If we use the results of previous subsections, we can write a more explicit formula 
for Mho using for \{mt) in ( pO]) the expression 



\{mt) 



+ 



3t 



LR 



[1 + 27g)tR] + 52\h - PZ^lr 
^^"''SiXHu{ML)tLR-f3?tn 



167r2 
3 



K 



167r2 

+ 2 

' 12XnuX'Hu)Ml ^^^^ 



tlR 



+ 



Ll67r2 



167r2 



tlR + 



with 



, Ml 



mt 



2 ' 



1 Ml 



, Ml 
In— ^ 

MI 



(32) 



(33) 



Once again, let us remark that in this formula the couplings are running MS parameters 
evaluated at the indicated scales (such scales are not explicit when different choices 
amount to three- loop effects). In particular, the mass parameters are defined by Mr = 
Mr{Mr), Ml = Ml{Ml) and mt = mt{mt) (here we write Wit for the SM running top 
mass). This is important to compare | p!l| , pT[ with results in OS scheme. The connection 
to OS parameters is dealt with in Appendix C. 
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To have a simple expression we have kept gt, in eq. (p^), evaluated at different 
scales in different terms. They are related to gtifrit) by 



gt{ML) = gtimt) + (3^2tL - d.gfiML) , (34) 



with Pg2 given in (p^ and Sig^^Mi) in (|T6D. We have also separated explicitly the two- 



loop threshold correction to Xh as 62\h- This correction can be most easily extracted 
from the two-loop effective potential and is presented in Appendix B. 

3 One-loop 'improved' formula? 



In the case of a single supersymmetric threshold, Ms ^ mt, it was shown in 1 14,15 



(following similar ideas already presented in ||T2[) that two-loop logarithmic corrections 
to M^o could be absorbed in a one-loop expression with running parameters evaluated 
at judiciously chosen scales. The compact formula obtained allowed a simple approxi- 
mation for the Higgs mass. Could a similar 'improved' one-loop expression for M^o be 
found for the case with a hierarchical stop spectrum? 

The derivation of the compact 'improved' formula in [jl5| is best understood in the 
RG approach. It starts with the expression 



X{mt) = \{Ms)- (3^^\mt) + (3^^ in lUf - l-f^ In^ Hlf + ... (35) 



2 din Q"^ mf 



which is the version of (^ which applies to the degenerate case Mr = Ml = Ms- As 
usual, given A(mt), the Higgs mass is obtained by (pO|). The idea is to use the freedom 
in choosing the scales in i3^\Tnt) ln[M|(M5)/mj (m^)], in (pSj), and lf{mt), in (|30|), to 
absorb two-loop corrections. 

The two-loop leading-log term in 



1 c/M^^ . iMl 

is easy to absorb by noting that 

SO that, if we choose = rritMs, the last term in ( |57D cancels (|5^. This choice of 
scale was first advocated in |]T2| and later confirmed by |T^,|TB|. 
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If we also want to absorb the two-loop next-to-leading-log term in (|35|) 



we have to choose appropriately the scales and Q'^ in 



/3iV(Q.)lni^^f^, (39) 



,(1)^2,^ MUMs) 
mKQ't) 

(which appears as a one-loop correction in M^o) to cancel ([38|). That a successful choice 
exists at all results from the happy accidental relation between the RG functions of the 
Standard Model: 



(40) 



This gives (see [|^0) Qt^ = M^rrit and Qj, = rrit. Further two- loop next-to-leading-log 
corrections to Mho associated with one-loop threshold corrections [like those coming 
from Xh{Ml) = X'^^Ml) + 6iXh{Ml) + if the parameters were evaluated at rrit 
instead of Mi] are trivial to absorb (in the previous example, just by choosing the 
scale to be Ml instead of rrit). 

The success of this program rests on two pillars. One is the existence of a single 
threshold at Ms, which allows a simple absortion of two- loop leading-log corrections. 
The second is relation (0). Note, however, that there is nothing fundamental about 
this relation. In fact, it is spoiled if one includes electroweak gauge couplings. Also, 
there is no reason to expect that leading or next-to-leading corrections beyond two 
loops will be given correctly by the advocated choice of scales. In conclusion, the 
compact formula was simply a useful approximation for the two-loop result, with no 
pretension to being fundamental. 

In view of the above discussion, the derivation of a compact one-loop 'improved' 
approximation to M^o in the case of a hierarchical stop spectrum looks problematic. 
To begin with, eq. ( ^OD holds for the running of the Higgs quartic coupling between 



Mji and rrit but not between Ml and M/?; that is, does not satisfy (40). Moreover, 
there are further complications associated with threshold corrections at M^j which were 
absent in the case of degenerate stop masses. All this implies that one cannot absorb 
two-loop next-to-leading-log corrections in any simple way. 



^We have left out of the discussion the comphcations associated to the fact that /jj^^^ contains a 
piece —XjHi which is formally of higher order. The effects of including this term properly can be 
reabsorbed in the scales at which one evaluates U^. The final result is as presented in |15 
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For two-loop leading-log corrections there are no threshold complications but the 
difficulties with two regimes of running (above or below M/j) remain. It is certainly 
possible to absorb these corrections using a trick similar to that in (^) in two separate 
one-loop terms 

-;3l"(Q,)ln^-<W.)lng, (41) 

with = Mimt and Q2 = Murrit, but /?Ah is a function of couphngs like Xhu which 
are not present in the low-energy theory and one would like to do better than that. An 
attempt in that direction was made in [Q, which for this type of spectrum advocates 
the use of 

-/3i"W3)ln^, (42) 

with the scale Q3 defined byQ 

|„4lnM£^?«=b.^££ + b.^£«, (43) 

mf mf mt mt 

as the one-loop redefinition of scales which absorbs two-loop leading-log corrections. 
In view of eq. (^Tj), for this approximation to be successful, (3^^ has to be related 
to (3^^ somehow. If we substitute Xhu in Px^ by its threshold value {^^, neglecting 
Xt-dependent terms (these could be absorbed eventually in Xj-dependent one-loop 
corrections) one has 

/J2 - itf ' , (44) 
which only holds a.t Q = M^. Were this to imply 

^ - 1^ (45) 

then, a one-loop approximation of the form (|^) could be devised, and in fact we would 
find a scale Q3 given precisely by (^3[) . However, (|4^) does not imply (|45|) , which in 
fact does not hold, as can be checked from the RGBs given in the previous subsections. 

The best we can do to absorb two-loop leading-log corrections in a one-loop formula 
is to use the freedom in choosing the scales of both mt [or, equivalently (3^^ 
although this only works for the Xj-independent corrections. If we focus only in the 
Xt = case, the one-loop 'improved' formula we find is 



4:7!-^v [Qv) mf 



^With some refinements like Ml — > M^^ and Mj^ M^^ which, in our approximations, do not 
affect the leading-log corrections we are discussing here. 
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with the scales Qt and Qy given by 

In — In — = - 

mf mf 6 



rrit rrit rrit rrit 

which differs from the Q3 defined in (|43|) , and 



(47) 



mt mf Mr 



Formula ( ^B]) successfully reproduces the two-loop leading logarithmic corrections to 
Mlo for Xt = 0. 

However, for the general case, with non-zero Xt, we conclude that there is no simple 
way of absorbing two-loop leading and next-to-leading logarithmic corrections to M^o in 
the case of a hierarchical stop spectrum. This does not preclude the possible existence 
of scale choices which minimize two-loop corrections, but we do not try to find them 
in this paper. 

4 Comparison with previous literature 

Several papers have studied two-loop radiative corrections for non-degenerate stop 
masses before. In the previous discussion, we have already mentioned ref. [jl2[, which 



performed a thorough study of radiative corrections to M^o using RG techniques to 
derive analytical and one-loop 'improved' approximations that take into account two- 
loop leading logarithmic corrections to the Higgs boson mass. As we have shown. 



our results do not agree with those of in the case in which we focus, with widely 
different stop masses. A possible reason for this discrepancy has been advanced above. 
Similar studies were performed also in ref. |]TT| , that uses a combination of effective 



potential and RG techniques to obtain the two-loop leading-log corrections to Mho 
working in the more general case of arbitrary pseudoscalar mass M^o (this complicates 
the analysis because, for low M40 , the low energy effective theory is a two Higgs doublet 
model). However, if we compare with the results of that paper for the case M^o = Ms 
we again find they disagree with our results, and the reason is similar to the one already 
mentioned: although at one-loop there is a simple relation between some RG-functions 
of the effective and full theories, its derivatives with the scale (which are necessary to 
get correctly the two-loop-leading-log terms) are more complicated than the one-loop 
relations suggest. 
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Before proceeding with the comparison to other previous analyses, we can aheady 
draw some imphcations of our resuhs. It is simple to show, focusing in two-loop- 
leading-log corrections and zero stop mixing for simplicity, that our M^o is higher than 
the previous estimates just commented ITHJl^ by the amount 



This formula assumes that the one- loop result is expressed in terms of rnt{'mt) and 
v{mt). As expected, the discrepancy (|49|) is larger the higher Ml/Mj^ is, disappears if 
Ml ~ Mr and is always positive. Numerically, it increases M^o up to ~ 5 GeV in the 
most extreme cases (say M^ — 2 TeV and M/j ^ m^). The first obvious implication 
of this upward shift of M^o is for experimental analyses of SUSY Higgs searches that 
made use of those previous mass calculations, whenever they were applied to scenarios 
with widely split stop masses. The same implications will also follow for theoretical 
studies in similar circumstances. 

Some two-loop radiative corrections to M^o (those which depend on the QCD gauge 
coupling) have been computed also diagrammatically [0, so we could make a partial 
check of our results. However, a complete expression for the diagrammatic result, 
applicable when the diagonal stop masses (M^ and Mr) are different, is too lengthy 
and has not been published. It would be interesting to make this comparison]]. 

Finally there are two-loop calculations of M/^o based on the use of the MSSM effec- 
tive potential. The first of them was the work by Hempfiing and Hoang in 0, which 
computed this potential for sin/5 = 1 and zero stop-mixing, extracting from it the 
radiative corrections to the Higgs mass. This work was extended later on by Zhang in 



13|, which added the QCD two-loop corrections to the effective potential for generic 



values of tan f3 and non-zero stop mixing. Finally, ref. [|T5| included also the top- Yukawa 
two-loop contributions to the MSSM effective potential. This last calculation agrees 
with previous effective potentials in the different limits in which those apply, so that 
we will only discuss here the comparison of our RG two-loop corrections with those 
that can be obtained from the effective potential as presented in ||15|| . 

The procedure used to get the Higgs mass from the MSSM two-loop effective poten- 



tial is similar to the one used and explained in [0,^. We first expand this potential, 

^The agreement between these two- loop diagrammatic results and the corresponding RG and/or 
effective potential results was first shown in . Unfortunately, this comparison had to be limited to 
the case of heavy and degenerate diagonal stop masses, AI^ = Mj^. = Msusy 3> m^. (See also p7[|) 
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V (which, for large M^o , is a function of the hght Higgs field through its dependence on 
rrit and the stop masses) in powers of mf/Ml, mf/M'^ and Mj^/Ml. Next we compute 

which gives an approximation to M^o as the second derivative of the potential (with 
respect to the Higgs field) at its minimum. Then, to get the physical Higgs pole mass, 
we correct for non-zero external momentum effects by adding to (^Dp the quantity 

H;,;,(0) - Uj^hiMlo) , (51) 

where Uhhip^) is the Higgs boson self-energy for external momentum p. What results 
from this procedure is an expression for M^o in terms of scale dependent DR '-running 
parameters. To make contact with our results in this paper we have to convert the 
parameters that survive below the SUSY threshold to the MS scheme, taking into 
account one-loop threshold corrections. For example, mt{Q), Mr{Q) and v{Q) are 
expressed in terms of its low-energy counterparts mt{Q), Mr{Q) and v{Q) (these 
relations are given in Appendix B). The final step is to express all running parameters 
in terms of their values at the scales at which we evaluate them in our RG approach: 
v{mt), mt = mtimt), Mr = Mr{Mr), Ml = Ml{Ml) and Xt(Mi). After doing this, 
we find an expression for M^o which is manifestly independent of the renormalization 
scale Q and which agrees exactly with our effective theory result (we give its expression 
explicitly in Appendix A). As a bonus, the effective potential calculation gives also the 
two-loop non-logarithmic correction, presented in Appendix B. 

This agreement is the best check of our results. It speaks greatly of the power of 
effective potential techniques: the effective potential V has built-in all the structure 
of RG-functions and threshold corrections that we had to compute afresh in the RG 
approach. Nevertheless our effort was not wasted because this structure which remains 
buried in V , needs to be made explicit to implement the resummation of logarithmic 
corrections to all loops of the RG programme. When the hierarchy in the stop masses 
is only moderate there is no need to resum large logarithms and one can revert to 
the use of the plain effective potential, which still gives correctly two-loop radiative 
contributions to M/jO. In other scenarios that may introduce large logarithms {e.g when 
stop mixing effects are large and cause a significant splitting of the stop masses or for 
more complicated patterns of SUSY particle masses) one should work out the relevant 
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effective theories and the corresponding RG-functions and threshold corrections. Still, 
to two-loop order the results of such calculations must agree with those coming from 
the effective potential in the same regime of parameters. 

Needless to say, it would be extremely useful to have the tools necessary to dig up 
all this structure from the effective potential directly. Studies on multi-scale effective 
potentials ^6], and in particular the recent proposal in [^, are a first promising 



step in that direction and it would be interesting to continue the development of such 
techniques. 

5 Conclusions and outlook 

The mass M^o of the MSSM light Higgs boson receives sizeable radiative corrections, 
the most important of which depend on the details of the stop spectrum (masses and 
mixing). In this paper we have revisited the calculation of the radiative corrections to 
MhO in the case of a hierarchical stop spectrum, M^^ ^ M^^ ^ with moderate stop 
mixing. We have used an effective theory approach to identify in these corrections non- 
logarithmic contributions (which can be interpreted as threshold corrections at different 
energy scales) and logarithmic contributions (which arise from renormalization group 
running of parameters between different energy thresholds). We have performed this 
calculation neglecting in radiative corrections all couplings other than the top Yukawa 
coupling and the QCD strong gauge coupling. Within this approximation we collected 
all radiative corrections to M^o up to, and including, two-loop terms. Our results 
correct previous calculations of two-loop leading-log corrections to the Higgs mass 



that appeared in the literature and are widely used |11,12], while we find complete 



agreement with the results of previous analyses based on effective potential techniques 



H| , p!5| . Numerically, we find that two-loop leading-log corrections increase M/jO by up 
to 5 GeV (in some cases with a large hierarchy of stop masses) relative to the results 
computed in |]TT],|T2]. This has obvious importance for the theoretical input used in 
experimental analyses. 

Beyond this comparison to previous findings, the results obtained can be used as the 
starting point of a numerical evaluation of the Higgs mass which can resum leading and 
next-to-leading logarithmic terms to all loops using renormalization group techniques. 
This resummation is mandatory to get an accurate calculation of M^o if the hierarchy 
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of masses in the stop sector is large. We defer such numerical analyses to some future 
paper. 

The present analysis can be extended in several ways. It is simple to add radia- 
tive corrections to Myfi that depend on the bottom Yukawa coupling, which may be 
important for large values of tan/5. Also, formulas similar to the ones derived in this 
paper could be found for the case of a light stop, with M^^ ~ ""^t- Other extensions 
of our results include the study of the radiative corrections to M^q when the hierachy 
between stop masses is due to a very large value of the stop mixing parameter, Xt or 
the connection of effective theory methods with the methods of multi-scale potentials 



developed in p6| , p7 |. 
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Appendix A: Explicit Higgs mass formula 



In this Appendix we write down the Higgs boson mass computed in the main text, 
exphcitly separating the radiative corrections as 

M^o = M| cos^ 2(3 + AillM^o + AinllM^o + A2llM^o + AajvLiM^o + A2nnllMIo , 

(A.l) 

where AillM^o contains the one-loop leading-log corrections, AinllM^o the one- 
loop next-to-leading-log corrections, A2llM^o the two-loop leading-log contributions, 
A2NLLMI0 the two-loop next-to-leading-log contributions, and, finally, A2nnllMIo, 
the next-to-next-to-leading-log (that is, the non-logarithmic) corrections. For two- 
loop terms, this division depends on the scales at which the parameters in one-loop 
contributions are evaluated. We take rrit = rnt{mt), v = Vt{mt), Mr = Mr^Mr), 
Ml = Ml{Ml), Xt = Xt{ML) and Yt = Yt{ML). Here mt(Q), v{Q) and Mr{Q) are 
MS -running parameters in the effective theories below Msusy, and Ml{Q), Xt{Q), 
Yt{Q) are DR '-running parameters in the full MSSM. For simplicity we take the stop 
mixing parameters to be real in the following expressions. 
For the one-loop pieces we obtain the well known results 
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and 



^nllMi^o 
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(A.3) 



where we have kept up to C(M|/M2) terms and used Xt = Xt/ML, % = Yt/ML. 
For two-loop corrections we find 
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with as = at^sM = 5't/(47r); and 
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with at = hf/{ATT). The two-loop non- logarithmic corrections {^2Nnll 
puted in the next Appendix. 



M^o) are com- 



Appendix B: Two-loop threshold corrections to the 
Higgs quartic coupling 

The threshold corrections to the Higgs quartic self-coupling, Ajy, appear in the match- 
ing of different effective theories at a given energy scale. They cannot be calculated 
by renormalization group methods but rather must be either computed by direct dia- 
grammatic calculations or extracted indirectly from the effective potential. The latter 
is the simplest method and is the one we follow in this Appendix. 

To illustrate the procedure, we consider first the simpler case of a unique super- 
symmetric threshold at Ms mt) which corresponds to the common mass scale of 
all SUSY particles, and obtain the threshold corrections to \h (both at Ms and mt) 
at two- loop order. As always in this paper we keep only radiative corrections which 
depend on the top Yukawa coupling, ht, and/or the strong gauge coupling constant, 
Qs- Doing this we will reproduce results already presented in |jl5[. We discuss later 
on the case of a hierarchical stop spectrum. For simplicity we take the stop mixing 
parameters to be real in this Appendix. 
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B.l Degenerate stops 

At the supersymmetric scale, M5, we match the MSSM theory (the 'effective theory' 
vahd above Ms) to the SM (the effective theory vahd below Ms). In the effective 
potential formalism we therefore have to match the MSSM potential [VMssAii^t)] to 
the SM potential [VsM{mt)]. As indicated, the Higgs field dependence of these functions 
appears always through the top quark mass, either directly or through the stop masses, 
which are now given by mf^ = M| + ± rritXt. The top quark mass is a scale- 
dependent parameter: is the MSSM DR '-running mass, while mt is the SM running 
mass in MS scheme. Both effective potentials [VMssMi^t) and VsMijnt)] are known as 
a perturbative expansion up to two-loops: V = V^^^ + V^^'^ + V^'^\ The expression for 
Vmssm can be found in [|I3],|T3|. The two-loop SM effective potential was first computed 
in It can also be extracted from Vmssm by keeping only the contribution of non- 
sup ersymmetric particles, trading nit by frit and adding an extra piece to account for 
the difference between DR and MS renormalization schemes. This extra piece reads 

{Wnr AV^l = -16g!mt{l - In ^) ■ (B.l) 

As we are interested in the quartic Higgs coupling only, we expand Vmssm — Vsm in 
powers of rrit/Ms (or rrit/Ms) and keep only the term that goes like the fourth power 
of the top quark mass. To proceed with the computation, we next convert nit to 
using 



ml{Q)=ml{Q)\l- ^' 



2 r 



s 



1 I . M^s MgXt 

Q' Ml 



3/i- 



+ — ^ 
327r2 



;i-?)(i-.nf)-i 



(B.2) 

This shift is a loop effect and therefore has different impact on the different contri- 
butions to the m^-term depending on their loop-order: first, it does not matter for 
the contributions coming from V^'^^ because the correction would be of three-loop or- 
der; second, the quartic couplings in VmIsm gauge couplings, and so, transforming 
MSSM parameters there to SM ones gives a loop correction which involves gauge cou- 
plings and we are neglecting such effects; finally for I^^^^-contributions the correction 
is a two-loop effect of the same order of the terms coming from V^'^'^ and must be kept. 
After this shift is performed, one simply extracts the threshold correction 62X from the 
term quartic in m^, with logarithms disregarded (they are renormalization effects) to 
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obtain 



(B.3) 



However this is not yet the final result quoted in |jT5|, eqs. (29-30) in that paper. 
In addition to the high-energy threshold correction at Ms there is a finite {i.e. non- 
logarithmic) low-energy correction coming from V^^]^ through 

2 



^th 



Am; 



d 



dmf 



2 ; ^SM ■ 



(B.4) 



[This kind of correction is zero for the QCD part of V^J^ or for 1//^^ 



SMi- 



The final result is the sum of the high energy threshold correction As plus the 
low-energy correction At from ( |B.4| ): 



A^^) M' 



A5 + At 



(B.5) 



with 



A, 



+ 



y^3y2 



-2Xt - 



7 ^ 1 ^ 

3 * 12 * 



6 *. 



16X,^ + fx 



+ 4 

25 



13 

y 



(3 + lQK){AXtYt + 1;') + (15 - 2AK)Xl 



/14 



— + (1 + AK)X^Yt + — + 2AK]XfY^' - — + m]XfY^ 



and, from ( p.4| ), 



A, 



V 3 



19 



12 



TT 



(B.6) 



(B.7) 



We have used = g'^/^An), at = h^/^Air), Xt = Xt/Ms, Yt = Yt/Ms and the constant 



K ~ —0.1953256. This result now reproduces exactly that presented in [15 



B.2 Hierarchical stops 

Now we first match at Ms the MSSM as high-energy theory to the SM + Ir as the 
intermediate-energy theory (valid below Ms = Ml). The procedure is similar to the 
one used in the previous case but with an expansion of VMssAiifnt) — VsM+hi^t) 
powers of rrit/ML, Wit/ML and M^i/Ml. Again we keep only the quartic power of the 
top quark mass in this expansion. Next, notice the argument rfit for the potential in 
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the intermediate energy regime: this is because the top quark mass is insensitive to 
the tn threshold and runs the same above and below it, that is, with SM RGEs. The 
shift from mj to mt in the one-loop potential contribution is now given by 



mi{Q)=mi{Q){l- 



9l n 



67r2 



+ ln 



Ml 



-2- 



M, 



G 



Mr 



2 In 



Ml 



-1 



We also have to convert from Mr{Q) above M^ to Mii{Q) below Ml: 



167r2 



127r2 



AM? ( In 



M? 



+M] 



24(M| - X^) [In 

2 „2 , l V>2 2 



R 



, Mf 



and from v{Q) above Ml to v{Q) below Ml'- 



v{Q) = v{Q) 



1 - 



647r2 



(B. 



(B.9) 



(B.IO) 



For simplicity, in the rest of the formulas wc remove the line from M/j in the under- 
standing that it is the parameter of the intermediate-energy theory. After throwing 
away logarithmic terms, one is left with the following threshold contribution (at the 
scale Ms) to the Higgs mass: 

r2 /A^2 X Q /o 9 . ,2 



A, 



+ 4Xf 



1 Y + y + (6 + TT )Xj + - 95JXj + (2 + TT )X^ Sf^ 



+ c 



3 2 



"3 



+ ^]X' 



+ — + 72KjXtYt 



+ -y/ + (62 - 7r^)Xf - 8(7 + 27K)X^Yt 
V9 3 ;** Vl8 3 



108K]X^Y^^ 



(B.ll) 



For the matching at the intermediate scale M/j we proceed in the same way with 
the expansion of V^M+tjil^t) ~ ^M(^t)- This expansion gives a threshold correction 



to M^o-. 



27r3^;2' 



87r3^;2 



:b.i2) 
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Finally, there is a low-energy correction identical to the one written in (|B.7| ). The 
end result for the two-loop finite correction to the Higgs mass is the sum of the three 
contributions just described: 



^th ^^-^ho 



A 



2NNLLMf^o 



As + A 



R 



At 



:b.i3) 



Appendix C: On-shell parameters 

We write down in the following the useful relations between the scale-dependent pa- 
rameters of the top-stop sector and the corresponding physical on-shell quantities, par- 
ticularizing general known results to our case. For simplicity we take the stop mixing 
parameters to be real in this Appendix. 
The tree-level mass matrix for stops is 



MhiQ) MUQ) 



Ml + mtXt 
mtXt Ml + ml 



{C.l] 



where we have neglected gauge couplings, which contribute to this matrix through 
D-terms, and indicated the implicit scale dependence of the matrix entries. The one- 
loop self-energy corrections for top and stops in the MSSM can be found in In 
the definition of physical parameters, care should be given to the choice of external 
momentum in these self-energies: the scale independence of the quantity in question 
crucially depends on that choice. For example, consider the momentum-dependent 
loop correction to the stop mass matrix (|C.1|) : 



AM 



(C.2) 



which also depends on the scale explicitly. Calling M? the radiatively corrected stop 
mass matrix, M? -|- AM?, we find the following scale dependence (at one loop): 

d 



167r 



1677^ 



-Tr M? 



hi{Mi^ + 2Ml^-2p' 



d 



rDet M? 



3Det m}-p\Ml^ + 2Ml^) 



(C.3) 



cilnQ2 

Using these equations, it is straigthforward to show that the momentum-dependent 



eigenvalues of M? being solutions of the equation 



x-xTi M| + Det M| 



(C.4) 
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are scale invariant if the external momentum squared is equal to the eigenvalue itself 
[i.e. on-shell). 

A stop mixing-angle {6t) which includes one-loop radiative corrections {6t) can be 
defined as the angle of the basis rotation which diagonalizes M? + AM?. As such, it 
depends on the value of the external momentum. The choice of that external momen- 
tum required to obtain a scale independent definition of the radiatively corrected stop 
mixing angle, 9t, is univocally fixed by demanding 



167r^ 



d 



dlnQ- 



tan 2^. 



, 2 (Mir + Mir. - 2p2) 

— ; tan2^t = , 



(C.5) 



which is satisfied for = {Tnf^ + m|^)/2. Once a scale-independent mixing angle 6t 
has been obtained in this way, we can define a physically meaningful, 'on-shell' mixing 
X^^, by the relation 

2MtXf^ 



sin26't 



12 



(C.6) 



Following these prescriptions, the final results for the physical parameters in the 
stop sector are as follows. For the on-shell stop masses we get 
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while for the physical stop mixing parameter we find 



X^' = Xt{Q) + 
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with 



^.^arctan^. 



(C.IO) 



Finally, the on-shell top-quark mass is 
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(c.ii: 



We also need to relate the running vev, f (Q), to some observable, like the mass of 
the W^. With the type of SUSY spectrum we consider in this paper, we find 
4 
7 



v\Q) 



^[M^ + Re n^^^(M^)] 



(C.12) 
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This vev corresponds to the minimum of the one-loop effective potential. With this 
definition (others are possible, see [Q) there are no explicit tadpole contributions in 
(|C.12D . A similar relation, but for the SM vev, v, reads 



4M^ 

g2 



3hfs 



327r2 
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(C.13) 



In addition, we must consider the relation between the Higgs pole mass, and the 
mass obtained from the effective potential. They are related by the shift 



-T{hh{Ml,) + T{hh{Q) 



^-^Mlo 
327r2 



Xt 



' , Ml 



Ml 
Ml 



(C.14) 

Note the partial cancellation that occurs when this correction is considered together 
with (^). 
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Some of the threshold corrections and RGEs we present in the main text can be 
obtained from the general expressions for self-energies in and RGEs in |^ and 



decoupling SUSY particles in them. (This works directly for two-point Green functions 
and through low-energy theorems for n-point ones). We have checked that, whenever 
applicable, our results agree with such alternative procedures. 
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